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ABSTRACT 


In  1964,  Schoeffler  [15]  developed  the  theory  of  continuous 
equivalent  network  by  applying  a  continuous  transformation  matrix  to  the 
nodal  admittance  matrix  of  a  given  passive  network.  He  applied  the  concept 
to  minimize  the  sensitivity  of  a  network  function  to  variations  in  the 
element  values  [16] .  More  recently  other  papers  on  this  subject  have 
appeared  [2,9].  All  these  authors  have  used  the  steepest  descent  method 
to  generate  the  equivalent  networks.  But  the  detail  of  the  computational 
procedure  is  not  described. 

In  this  thesis  the  concept  of  continuous  equivalence  is  used  to 
optimize  networks  with  respect  to  two  performance  criteria,  namely,  distri¬ 
bution  of  elements  and  sensitivity.  The  transformation  matrix  is  applied 
to  the  loop  impedance  matrix  instead  of  the  nodal  admittance  matrix  as  done 
by  Schoeffler  in  his  original  work.  One  of  the  features  of  the  method  of 
solution  adopted  by  Schoeffler  and  others  is  the  need  to  assign  new  values 
for  the  elements  of  the  transformation  each  time  the  generating  equations 
are  solved.  Since  these  element  values  are  arbitrary  this  task  could  be 
tedious.  Although  the  steepest  descent  method  is  also  used  in  this  thesis, 
a  way  has  been  found  by  means  of  a  substitution,  which  enables  the  computer 
to  automatically  assign  the  appropriate  new  values  to  the  elements  of  the 
transformation  matrix  each  time  the  generating  equations  are  solved.  Numeri¬ 
cal  examples  of  one-port  and  two-port  networks  are  solved.  Certain  limita¬ 
tions  of  Schoeffler's  continuous  equivalence  concept  are  discussed  at  the  end 
of  this  thesis. 
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Chapter  1 


EQUIVALENT  NETWORKS 

1-1  Definition :  Two  networks  N^  and  N£  are  said  to  be  equivalent 

if  the  driving  point  and/or  transfer  immittance  of  N^  is  the  same  as  that 
of  N^.  The  values  of  some  or  all  the  elements  of  N^  may  be  different  from 
those  of  N£.  This  definition  will  be  used  in  this  thesis.  Several  other 
connotations  can  be  given  to  the  word  "equivalence."  We  will  also  be 
concerned  only  with  passive  networks.  Fig.  1-1  represents  two  networks, 
whose  driving  point  impedances  are  Z^Cs)  and  Z^C3)  respectively,  s  is  the 
Laplace  Transform  variable. 


are  the  elements  connected 
inside  the  box. 


Fig.  1-1  Two  one-port  passive  networks. 
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These  two  networks  will  be  equivalent  if 

Z  (s)  =  Z,  (s). 
a  b 


1-2  Howitt's  Approach. 

T 

Howitt  [5]  propounded  the  theory  of  equivalent  networks  in  1931. 
But  two  years  before  him  Cauer  [4]  had  already  propounded  the  same  theory. 
Their  philosophy  is  based  on  a  long  known  "fact"  but  to  which  nobody  before 
them  paid  any  serious  attention.  The  "fact"  is,  for  a  given  network  with 
specified  values  of  elements,  one  and  only  one  impedance  function  can  be 
realized;  but  for  a  given  impedance  function,  there  may  be  an  infinite 
number  of  networks  that  can  be  generated. 

•j—j* 

Howitt  uses  a  form  of  congruent  transformation  matrix  A,  on  the 
loop  impedance  matrix  of  a  given  network  thereby  generating  a  family  of 
networks.  The  elements  of  the  transformation  matrix  A  are  picked  arbi¬ 
trarily  so  that  all  the  new  networks  have  the  same  driving  point  impedance 
as  the  given  network. 

The  transformation  matrix  A,  will  be  referred  to  as  Cauer 
transformation  matrix  in  this  thesis,  and  it  is  of  the  form 


A  = 


(1-1) 


Numbers  in  square  brackets  refer  to  articles  or  books  listed  under 
Bibliography . 

TTSymbols  with  a  bar  underneath  refer  to  vectors  or  matrices  as  the  case 
may  be. 


. 


> 
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in  order  to. keep  the  driving  point  impedance  at  certain  ports  constant. 
_u  is  the  identity  or  unit  matrix. 


The  choice  of  A  rests  on  a  theorem  which  follows: 


VS> 


Theorem.  Two  (n  +  m)  x  (n  +  m)  loop  impedance  matrices 
and  Z^(s)  have  the  same  n  x  n  terminal-impedance  matrix  Z,(s) 

Zf  =  &Z.  A  . 


if 

(1-2) 


with  A  = 


|\  ! 
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and  &  is  non-singular. 

Note  that  J\,  can  be  made  equal  to  A*";  the  transpose  of  the 

matrix  A. 

In  his  work,  Howitt  makes  A  a  constant  matrix.  And  by  choosing 
—1*  —2  some  ar^itrary  fashion  he  arrives  at  a  family  of  networks  some 


* 

Proof  of  this  theorem  is  in  Appendix  I. 


- 


. 
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of  which  have  the  same  number  of  elements  as  the  given  network  and  some  of 
which  have  fewer  elements. 

1-3  Schoeffler *s  Theory  of  Continuously 

Equivalent  Networks 

Recently,  Schoeffler  [15]  has  improved  on  Howitt’s  theory  by 
making  the  transformation  matrix  a  function  of  some  dummy  variable  x  and 
then  applying  the  transformation  directly  to  the  elements  of  the  network. 
This  results  in  the  elements  of  the  network  being  functions  of  x.  At  each 
value  of  x,  the  new  network  is  equivalent  to  the  original  or  given  network. 

Schoeffler  developed  his  theory  for  the  admittance  case  (i.e., 
all  elements  are  connected  in  parallel) .  A  summary  of  this  theory  is  given 
in  Appendix  II.  In  practice,  networks  with  elements  connected  in  series 
occur  frequently.  Direct  application  of  Schoeffler 's  results  is  not  found 
to  be  very  convenient  in  view  of  the  large  number  of  equations  which  may 
have  to  be  solved.  Hence  it  is  found  desirable  to  derive  parallel  results 
for  the  impedance  case;  i.e.,  the  case  in  which  the  elements  of  the  network 
are  connected  in  series.  This  derivation,  which  will  be  presented  below, 
follows  Schoeffler’s  development  very  closely;  the  difference  occurs  in 
the  type  of  restrictions  placed  on  the  transformation  matrix,  TT. 

Suppose  we  are  given  a  network  which  can  be  described  by  the 

equation : 

V 

where  V 

I 


=  Z_  I  .  (1-3) 

-  tV  V2 . V*  .  (1-4) 


- 

. 


' 


' 


’ 
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and  z_  is  the  n  x  n  loop  or  mesh  impedance  matrix. 

The  branch  impedance  matrix  Z,  of  the  network  in  general  is 

— b 

of  the  form 


Vs)  -  h  +  ah  +  i£ 


(1-6) 


where 


4 


is  the  branch  resistance  matrix 


L,  is  the  branch  inductance  matrix 
— b 

Cb1  is  the  branch  inverse  capacitance  matrix 


If  we  assume  that  there  is  no  coupling  between  any  two  elements,  R^, 


-1 


and  are  all  diagonal  matrices.  For  example: 


R1  0 . 0 

o  r2 . 9 

R3 

0  0 . -  R 


where  m  is  the  number  of  branches  in  the  network. 

The  loop  impedance  matrix  Z_  can  be  expressed  in  terms  of  the 
branch  impedance  matrix  by  the  equation: 


Z  =  Bt  Z,  B 

—  —  — D  — 


(1-7) 


where  B  is  the  branch-mesh  incidence  matrix  and  the  elements  of  15  are 


defined  by 


f  1  If  branch  j  is  traversed  by  mesh  k  in 
the  branch  direction. 

b,,  =  ~  -1  If  branch  j  is  traversed  by  mesh  k  in 

J  the  opposite  direction. 

\^0  If  branch  j  is  not  traversed  by  mesh  k, 


Using  Cauer  transformation  A,  on  the  loop  impedance  matrix 
we  get  another  loop-impedance  matrix  Z*  given  by: 


Z'  =  A  Z  A 


(1-8) 


Substituting  equation  (1-7)  into  (1-8)  we  get 


t  „t 


Z'  -  A  B  Z_  B  A 

—  —  —  — D - 


(1-9) 


If  a  transformation  is  now  applied  to  the  branch  impedance 
matrix  Z^,  we  get  a  new  branch  impedance  matrix  _Z ' v  given  by 


Z'  =  T  Z  T 
—  b  —  — b  — 


(1-10) 


A  new  loop  impedance  matrix  7/  can  also  be  formed  from  equation 
(1-10)  using  equation  (1-7) . 


Thus : 


Z’  =  BC  TZ 


B  T  I  T  B 


(1-11) 


For  Z'  in  equations  (1-9)  and  (1-11)  to  be  equal 


. 


' 


■ 


- 


' 


- 


'  1 
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OR 


(1-12) 


B  A 


T  B 


Equation  (1-12)  describes  the  restriction  that  should  be  placed  on  T?  if 
an  equivalent  network  is  to  result. 

Since  the  branch  impedance  matrix  Z,  can  be  separated  into  the 
different  element  kinds  as  represented  by  equation  (1-6) ,  it  follows  that 
the  transformation  matrix  can  be  applied  to  each  element  kind  separately 

without  affecting  the  overall  result. 

If  we  now  assume  that  the  transformation  matrix  T  is  a  function 
of  a  dummy  variable  x  the  resulting  elements  from  the  transformation  are 


then  functions  of  x,  i.e.. 


Let 


T 


u  +  AT 


(1-13) 


where  u  is  a  b  x  b  unit  matrix,  b  being  the  number 


of  branches 


Then 


T 1 1  Zfe  T 


— b  — b  —  —  — b 


Z,  +  Z,  AT  +  ATC  Z,  +  ATC  ZL  AT 


(1-14) 


If  we  let 


AT 


2l]_  Ax 


where  Ax  is  a  scalar  and 


T^  is  of  the  form  (for  a  four  branch  network) 


'  •• 

- 

IF 


Il 


then 


Now  let 


8  - 

C2 

C3 

C4 

C5 

C6 

*7 

t8 

C9 

t10 

Cn 

t12 

t 

13 

t 

14 

£ 

15 

t 

16 

T,  Ax  +  T 

3  — i  — 

th 

Ax  + 

*  Il 

(x)  , 

h 

-  h 

-1  =b  -1 


and 


z;  = 


Z^(x  +  Ax) 


(1-15) 


equation  (1-15)  becomes 


Z,  (x  +  Ax)  -  Z,  (x) 

— - 12 -  =  4U)  T^x)  +  .T*(x)  Z^x) 


Ax 


+  T£(x)  Z^(x)  Tx(x)  Ax 


(1-16) 


Taking  the  limit  in  equation  (1-16)  as  Ax  -►  0  we  have 


dA 

dx 


4  ii  +  ii  4 


(1-17) 


The  problem  now  is  to  choose  the  elements  of  T^  so  that 
equivalent  networks  will  result.  But  even  before  we  choose  the  elements, 
we  have  to  place  a  restriction  on  T^  which  is  similar  to  that  placed  on  T_. 


, 


' 


,1 


■ 


-  9  - 


Substituting  equation  (1-13)  for  T*"  in  equation  (1-12)  yields 

B*  =  BC  (u^  +  Ax)  .  (1-18) 

Cauer  Transformation  A  is  also  written  in  the  form 

A  =  u_  +  A  A  =  _u  +  Ax  .  (1-19) 

Substituting  (1-19)  into  (1-18)  we  have: 

A^  B*  =  BC  .  (1-20) 

Equation  (1-20)  gives  the  restriction  on  T^. 

A  few  remarks  on  the  nature  of  the  matrix  A^  may  be  appropriate 
here.  Suppose  that  it  is  required  to  keep  the  driving  point  impedance  at 
terminal  pairs  1-1’  and  2-2’  of  the  two-port  network  in  Fig.  1-2 
constant.  A  look  at  the  proof  of  the  theorem  in  Appendix  I  immediately 
shows  that  if  the  driving  point  impedances  at  certain  terminal  pairs  are 
kept  constant,  the  transfer  impedances  between  those  pairs  of  terminals  will 
also  be  constant. 


1  O 


l’o 


C  2 


Fig.  1-2  Two-Port  Passive  Network. 
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A  comparison  of  equations  (1-19)  and  (Al-3)  shows  that  for  the 
above  conditions  to  be  satisfied,  the  rows  of  A^  corresponding  to  those 
terminal  pairs  must  have  zero  elements. 

For  example,  if  there  are  four  loops  in  Fig.  1-2,  then  A^  will 
be  of  the  form 


(1-21) 


The  free  elements,  a^,  a a^  .  can  be  chosen  arbitrarily 

to  arrive  at  equivalent  networks.  However,  it  is  important  to  note  that 
these  elements  should  be  such  that  the  Cauer  transformation  matrix.  A, 
is  non-singular.  For  example  in  equation  (1-21)  if  a^  =  -1,  then  a^ 
must  be  non-zero,  to  ensure  that  A  is  non-singular. 

We  will  conclude  this  section  by  arriving  at  a  differential 
equation  for  each  kind  of  element  in  a  given  network.  Equation  (1-17)  was 
arrived  at  by  applying  the  transformation  to  the  branch  impedance  matrix 
.  If  we  follow  a  similar  procedure  and  apply  a  transformation  T^R,  T^l 
and  T.  to  R^,  and  respectively,  we  will  arrive  at  the  following 


equations . 


' 


■ 

i  i  ■  '-  ‘  fa  r 

■ 

* 

- 


11  - 


dx 


hhi 


+  hth 


+ 4  ^ 1 


We  also  obtain  a  relation  similar  to  equation  (1-20) : 


-1 


(1-22) 


(1-23) 


(1-24) 


(1-25) 


The  three  differential  equations  (1-22) ,  (1-23)  and  (1-24) ,  when 
solved  subject  to  equation  (1-25),  will  generate  the  element  values  as 
functions  of  x.  For  different  values  of  x,  a  family  of  equivalent  networks 
will  result.  Therefore,  these  differential  equations  will  be  referred  to, 
as  "generating  equations"  in  the  rest  of  this  thesis. 

If  any  of  the  element  values  should  become  zero,  during  the  course 
of  the  solution,  additional  constraints  should  be  imposed  on  the  generating 
equations  to  prevent  the  element  values  from  becoming  negative.  These  con¬ 
straints  will  also  hold  if  in  the  given  network  an  element  is  absent  (with 
respect  to  the  general  network)  and  we  want  to  keep  it  so.  The  nature  of 
the  additional  constraints  will  be  discussed  later  (Chapter  3,  example  3-1). 
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1_4  Application  of  Equivalent  Networks 

An  obvious  question  which  comes  to  one’s  mind  is,  Why  generate 
equivalent  networks?  The  need  for  generating  equivalent  networks  arises 
from  the  fact  that  circuit  designers  aspire  to  attain  certain  goals.  Some 
of  these  goals  include  circuits  with  minimum  number  of  elements,  with  best 
possible  element  distribution  and  with  minimum  sensitivity  of  the  immittance 
function  to  changes  in  element  values.  In  other  words,  circuit  designers 
like  to  ensure  that  the  circuits  they  have  designed  to  do  a  job  are  the 
"best"  ones. 

Hence,  instead  of  choosing  the  free  elements  of  A^  arbitrarily 
(which  amounts  to  choosing  T_  arbitrarily)  we  can  choose  them  to  attain  one 
or  combinations  of  the  above  goals.  Thus  the  main  application  of  the  con¬ 
cept  of  continuous  equivalence  is  in  the  optimization  of  networks.  However, 
it  can  also  be  used  to  eliminate  coupling  between  elements  as  well  as  to 
eliminate  negative  elements  from  a  given  network. 

1-5  Scope  of  This  Thesis 

In  this  thesis,  the  theory  of  continuously  equivalent  networks 
will  be  used  to  optimize  networks.  The  network  optimization  will  be  carried 
out  for  two  performance  criteria,  better  element  distribution  and  minimum 
sensitivity.  Examples  of  one-  and  two-port  networks  will  be  considered. 

The  limitations  of  continuous  equivalence  as  far  as  this  author  could  see 
will  also  be  discussed. 

One  of  the  contributions  of  this  thesis  is  in  the  procedure  used 
to  solve  the  differential  equations  which  result  from  the  application  of  the 
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equivalence  concept.  As  will  be  seen  in  the  next  chapter  these  equations 
contain  arbitrary  coefficients.  Schoeffler  seems  to  imply  that  values  have 
to  be  assigned  to  these  coefficients  repeatedly  to  satisfy  certain  equations 
thus  leading  to  solving  the  differential  equations  repeatedly.  A  more 
systematic  method  of  solving  the  differential  equations  is  presented  in  this 


thesis . 


W*S??V 


Chapter  2 

NETWORK  OPTIMIZATION— I  (Element  Distribution) 


2-1  Introduction:  As  was  pointed  out  in  the  previous  chapter,  the 

continuous  equivalence  concept  enables  us  to  generate  a  family  of  networks 
which  have  some  common  characteristic,  for  instance  the  same  driving  point 
impedance.  In  obtaining  these  new  networks  we  may  also  be  interested  in 
improving  or  optimizing  the  network  with  respect  to  some  specific  criterion. 
For  example,  consider  the  n-port  network.  Fig.  2-1. 


i 


1’  & 


Fig.  2-1  N-Port-Network 
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Let  us  assume  that  it  is  required  to  have  a  better  distribution 
of  values  in  all  the  elements.  As  in  all  optimum  problems,  this  will  have 
to  be  expressed  in  a  mathematical  form.  If  we  represent  the  criterion  by 
then  it  can  be  expressed  as 

MNP 

*  -  *R  y.  (R±  ■ R)2  +  \  f  u-i  -  v2  +  Ah  2l  (hi  ■  w2 

i=l  i-1  i»l 

.  (2-1) 

where  h.  =  i  and 
i  C. 

1 

AR,  A^  and  A^  are  the  weighting  factors  for  the  resistances, 

inductances  and  inverse  capacitances  respectively,  and  R,  L  and  h  are  the 

mean  values  of  the  R’s,  L's  and  h's  respectively.  For  instance 

R..  +  R«  + . + 

R  -  — - - - 

M 

Obviously,  $  is  a  measure  of  the  deviation  of  the  actual  element  values 
from  their  average.  Our  objective  will  be  to  minimize  $.  The  actual 
expression  chosen  for  $  is  only  one  of  several  ways  of  defining  the 
criterion. 

In  the  remaining  part  of  this  chapter,  two  examples  will  be 
solved  to  illustrate  the  optimization  with  respect  to  element  distribution. 
The  first  example  deals  with  a  one-port  network  while  the  second  deals  with 
a  two-port  network.  For  the  sake  of  simplicity,  we  will  attempt  to  get  a 
better  distribution  of  values  for  the  resistances  only.  This  means  that 


* 


' 


1 
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equation  (2-1)  simplifies  to: 

M 

$  =  (R.  -  R)2  .  (2-2) 

i=l 


The  examples  will  be  solved  in  two  stages.  First,  the  generating 
equations  for  continuous  equivalence  will  be  obtained.  Then  the  performance 
criterion  given  in  equation  (2-2)  will  be  minimized  using  the  generating 
equations  as  the  constraints.  The  method  of  steepest  descent  will  be  used 
for  arriving  at  the  optimal  solution. 


2-2 


Fig.  2-2  R-L-C-Network. 

Dimensions  are  in  ohm,  henry  and  inverse  farad  or  daraf. 


' 
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Consider  the  network  shown  in  Fig.  2-2.  The  values  of  the 
elements  are: 

Rx  =  3  R2  =  2  R3  =  4 

L1  =  1  L2  =  2  L3  =  3 

h1  =  2  h2  =  1  h3  =  0.5 

There  is  no  mutual  coupling  between  the  inductances. 

From  Fig.  2-2  we  can  write  the  branch  impedance  as: 


4 


R,  +  sL,  +  — 
1  Is 


0 


0 


0 


R0  +  sL0  H - 

2  2  s 


0 


0 

0 


R„  +  sL_  +  — 
3  3  s 


(2-3) 


— 

— 

R1 

0 

0 

L1 

0 

0 

i — 1 

X 

0 

0 

0 

R2 

0 

+  s 

0 

L2 

0 

+  i 

s 

0 

h2 

0 

0 

0 

R3 

0 

0 

L3 

0 

0 

h3 

— 

.  (2-4) 

The  branch-mesh  incidence  matrix  B,  by  inspection  of  Fig  (2-2) 


is, 


(2-5) 


0 


1 


i 


. 


The  generating  equations  are  the  same  as  equations  (1-22) ,  (1-23) 
and  (1-24).  For  convenience,  the  equations  are  written  below: 


d4 


dx 


...  [1-22] 


dx 


ib  In  +  IlL  h 


..  [1-23] 


dll 

-  4  lie  +  he  “b  •••  U-24] 


dx 


where  T^,  T^  and  _T  are  of  the  form 


1R 


— 

*1 

t2 

t3 

~i 

t2 

e3 

C4 

C6 

T  = 

1L 

C6 

C7 

t8 

C9 

*7 

t8 

_ 

— 

... 

T 

-1C 


t"  t"  t" 

1  z2  c3 

t"  t"  t" 

4  C5  6 

t"  t"  t" 

7  8  9 


It  must  be  remembered  that  for  equivalence  to  result  equation 
(1-25)  must  also  be  satisfied.  In  order  that  the  driving  point  impedance 
remains  invariant  under  the  transformations,  A-^  should  be  of  the  form 


I  J  k&Jm 

■  II  1  I  .1,,  J  I 

■ 


. 


, 

*■ 
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0  0 


Substituting  for  and  in  equation  (1-22)  we  have: 


R1 

0 

0 

2Rltl 

R,  t „+R„ t, 
12  2  4 

Rlt3+R3t7 

0 

R2 

0 

= 

Rlt2+R2t4 

2R2t5 

R2t6+R3t8 

0 

0 

R3 

R3t7+Rlt3 

R2t6+R3t8 

2R3t9 

— 

— < 

which  gives 


dR, 


dx 


■  2Rici 


dR, 


dx 


‘  2R2t5 


dR, 


dx 


2R3t9 


and 


R,  t  +  R0t.  =  0 

12  2  4 


V3  +  V?  “  0 


R2t6  +  R3t8  0 


Substituting  for  A^,  _B  and 


(2-6) 


(2-7) 


in  equation  (1-25)  and  substituting  the 


. 


' 

1  . 
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result  in  equation  (2-7)  the  following  relationships  are  obtained. 

2Rlt1  =  -  R2(a2  +  aP  +  R3ax 

2R2t5  =  R2(a2  -  ax)  -  R3ax  (2-8) 

2R3t9  =  R2(a2  +  a1)  +  R3(2a2  +  a^ 

If  equation  (2-8)  is  now  substituted  in  equation  (2-6)  the  result  is 


dR, 


dx 


=  -  R2(a2  +  a1)  +  R3a1 


dR, 


dx 


=  R2(a2  -  a^  -  R3ax 


(2-9) 


dR, 


dx 


=  R2(a2  +  a^)  +  R3(2a2  +  a^ 


Going  through  a  similar  procedure  for  the  inductance  and 


capacitance  group  we  have 


dL 


1 


dx 


=  -  L2 (a2  +  Zj)  +  R3a1 


dL, 


dx 


L2^S2  al^  L3al 


.  (2-10) 


dL, 


dx 


=  L2(a2  +  ax)  +  R3(2a2  +  a^ 


■ 


w*  ••  .  ' 

' 
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=  -h2(a2  +  a±)  +  h3ax 


h2(a2  a2_)  h3ai 


h2(a2  +  ax)  +  h3(2a2  +  a^ 


(2-11) 


In  equations  (2-9) ,  (2-10)  and  (2-11)  and  a2  are  not  known 
and  can  be  assigned  values  arbitrarily.  However,  in  this  example,  the 
values  will  be  chosen,  keeping  in  mind  that  the  distribution  of  values  of 
the  resistances  must  be  optimum. 

For  this  example,  equation  (2-2)  becomes 


$ 


. ..  [2-2] 


_  R,  +  R2  +  R3 

where  R  =  - 

3 

or  $  =  |(R^  +  R2  +  R3  -  R1R2  -  R1R3  -  R^)  .  (2-12) 

The  method  of  steepest  descent  will  be  used  to  determine  the 
values  of  x  which  ensure  that  $  goes  rapidly  to  its  minimum  or  at 
least  approaches  its  minimum  fast. 


d$ 

dx 


0  . 


For  $  to  approach  its  minimum  rapidly,  it  is  required  that 


■ 


- 


. 
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But  - —  = 


d$ 

dx 


3$_ 

3R-, 


dR 


dx 


1  3$ 


3R, 


dR_  dR„ 

_ 2  3$  .  3 

dx  3R_  dx 


(2-13) 


Evaluating 


3$ 

3R, 


3$ 

3R, 


and 


3$_ 

3R, 


by  differentiation  of  equation  (2-12) 


1  w“2  ”“3 
and  substituting  the  results  into  equation  (2-13)  and  then  substituting 


dR]_  dR2  dR3 

for  fa~  >  ~fa~  and  •  from  equations  (2-9),  (2-10)  and  (2-11)  into 
equation  (2-13) ,  we  get 


d$  2  o  ? 

fa  =  3  (a2(-4R1R2  -  2RxR3  +  2R2  +  4R^) 

+  a1(-2R1R2  +  2R1R3  -  2R2  +  2R3)}  .  (2-15) 

Let  us  define  a^  and  a2  by 

a2  -  -  -|(-4R1R2  -  2R1R3  +  2R2  +  4R3)  .  (2-16) 

a1  =  -  ■|(-2R1R2  +  2R][R3  -  2R2  +  2R3)  .  (2-17) 


Thus 


hence  is  always  negative. 

By  making  the  definitions  described  by  equations  (2-16)  and 
(2-17)  we  avoid  the  uncertainty  involved  in  assigning  values  to  the 
coefficients  a^  and  a2  and  repeatedly  solving  the  differential  equations 
until  we  arrive  at  the  optimum  values  of  a^  and  a2> 

For  this  choice  of  values  of  a^  and  a2  given  by  equations  (2-16) 
and  (2-17)  the  linear  generating  equations  (2-9),  (2-10)  and  (2-11)  become 


nonlinear  and  given  by 


' 

■ 

' 


(  >:i  f 
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d4 


dx 


“4 


S. 

dx 


M  L, 
- b 


% 

dx 


=  M 


4 


(2-18) 


(2-19) 


(2-20) 


where 


0  . 

|(6R^  -  6RxR2) 

|(2R1R2  -  2RxR3  +  2R2  -  2R3) 

M  = 

0 

|(2R1R2  +  4RxR3  -  4R2  - 

2R3> 

■|(2R1R3  +  2R3  -  2RxR2  -  2R2) 

0 

|(6R1R2  -  6R3) 

|(10R1R2  +  2R1R3  -  2R2  -  10R3) 

and  as  before 

4  “  [Ri  4  R3]t 

4  "  [L1  L2  L3]t 

4  =  [hl  h2  h3]C 


The  nonlinear  differential  equations  (2-18) ,  (2-19)  and  (2-20) 
are  solved  on  the  digital  computer  (IBM  360)  .  The  solutions  are  displayed 
in  Appendix  IV. 

Assuming  that  at  x  =  0.030  we  are  satisfied  with  the  element 
values,  Table  2-1  shows  the  comparison  between  the  new  element  values  and 
the  given  ones  as  well  as  the  respective  values  of  $ . 


. 


' 
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ORIGINAL 

NETWORK 

x  =  0 

NEW 

NETWORK 

x  =  0.030 

R1 

3 

2.921 

R2 

2 

2.790 

R3 

4 

2.860 

L1 

1 

1.056 

L2 

2 

2.664 

L3 

3 

2.005 

hl 

2 

2.151 

h2 

1 

1.186 

h3 

0.5 

0.215 

$ 

2 

0.008 

TABLE  2-1  Summary  of  Results 


Fig.  2-3  Variation  of  Resistive  Elements  with 
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x 


Fig.  2-4  Variation  of  Inductive  Elements  with 


OR 
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Fig.  2-5  Variation  of  Capacitive  Elements  with 
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If  we  plot  the  values  of  each  element  as  a  function  of  x,  the 
resulting  curves  may  be  used  if  necessary  to  obtain  approximate  analytical 
expressions  for  the  elements  as  functions  of  x  by  using  curve  fitting 
procedures.  In  Figs.  2-3,  2-4  and  2-5,  the  variations  of  R’s,  L’s  and  h's 
with  x  are  shown. 

Fig.  2-6  shows  the  new  network. 


2.921  1.056  2.151 


2.860  2.005 


0.215 


Fig.  2-6  Network  Equivalent  to  Fig.  2-1 


The  driving  point  impedance  of  the  original  and  the  new  net¬ 


works  are: 


The  expression  for  each  Z(s)  has  been  normalized  to  make  the  coefficient 
of  the  highest  power  of  s  equal  to  unity. 


' 
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Original  network  (Fig.  2-1) 


Z(s)  = 


>4  +  3.182  s3  +  3.770  s2  +  1.953  s  +  0.318 

s (0 .455  s2  +  0.540  s  +  0.136) 


New  network  (Fig.  2-6) 

n  s4  +  3.199  s3  +  3.690  s2  +  1.970  s  +  0.318 
Z(s)  =  - - - 

s (0.454  s  +  0.550  s  +  0.137) 

Table  2-2  shows  the  percentage  deviation  of  the  coefficients  of 
the  powers  of  s  for  the  Z(s)  of  the  new  network  from  the  corresponding 
coefficients  of  the  original  network. 


COEFFICIENTS 

%  DEV] 

[AT ION 

OF 

NUMERATOR 

DENOMINATOR 

4 

s 

0 

3 

s 

-  0.220 

0.250 

2 

s 

2.120 

-  0.732 

1 

s 

-  0.895 

-  0.514 

0 

s 

0 

TABLE  2-2  Percentage  Deviation  in  Coefficients 
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Example  2-2 

We  will  now  consider  the  optimization  of  a  2-port  network  using 
the  same  performance  criterion,  namely  element  distribution.  The  given 
network  is  shown  in  Fig.  2-7.  The  driving  point  impedances  at  both  ports 
will  be  made  invariant  under  the  transformation.  The  number  of  loops  being 
greater  than  the  number  of  ports,  is  by  no  means  an  accident.  That  this 
indeed  is  a  requirement  follows  from  the  proof  of  the  theorem  stated  in 
Chapter  1.  This  point  will  be  discussed  in  some  detail  in  Chapter  4. 

The  branch  impedance  matrix  Z^s)  an<*  t*ie  branch-mesh  incidence 
matrix  B  ,  are: 


Vs>  - 


r— i 

0 

0 

0 

0 

L1 

0 

0 

0 

0 

0 

R2 

0 

0 

0 

0 

L2 

0 

0 

0 

0 

0 

R3 

0 

0 

+  s 

0 

0 

L3 

0 

0 

0 

0 

0 

R4 

0 

0 

0 

0 

L4 

0 

0 

0 

0 

0 

R5 

0 

0 

0 

0 

L5 

___ 

(2-21) 


B  = 


1 

1 

0 

0 

0 


0 

0 

0 

1 

1 


0 

1 

-1 

-1 

0 


(2-22) 
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The  matrix  is  given  by 


The  generating  equations  are: 


(2-24) 


^4 

dx 


4  i1R  +  I JR  4 


(2-25) 


dx 


h  hi  +  hi  h 


(2-26) 


where 


and 


are  given  by 


— 

h 

t2 

C3 

C4 

C5 

C6 

C7 

t8 

t9 

C10 

T  = 

-1R 

C11 

t12 

t13 

t14 

t15 

C16 

t17 

t18 

t!9 

C20 

t21 

t22 

C23 

fc24 

C25 

”i 

t2 

t3 

*4 

‘f 

C6 

*7 

t8 

t9 

cio 

T  = 

-1L 

hi 

t12 

t13 

t14 

‘is 

t16 

fci7 

t18 

t19 

t20 

C21 

t ' 

22 

t23 

fc24 

C25 

— 

(2-27) 
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As  before,  equation  (1-25)  must  also  be  satisfied,  if  equivalent 
network  is  to  result. 

Thus:  Bfc  T?  =  A,  Bt  =  Bt  t5_  .  (2-28) 

Since  the  procedure  used  here  is  the  same  as  that  used  in  example 
2-1,  most  details  will  be  omitted.  The  important  equations  and  results 
will  be  simply  stated. 

Thus,  the  generating  equations  become 


dR, 


dx 


R2 (ai  -  a2  “  a3)  “  R3ax 


dR, 


dx 


R2(a1  +  a2  +  a3^  +  R3ai 


dR, 


dx 


R2(a3  -  a1)  +  R^ (2a^  +  a2  -  a.^)  +  R^(a3  +  a2) 


dR, 


dx 


-R„a0  -  R. (a,  +  a_  -  a0) 
3  2  4  1  2  3 


dRr 


dx 


R3a2  +  R4(ax  “  a2  “  a3> 


(2-29) 


dL. 


dx 


L2^al  a2  a3^  L3al 


dL, 


dx 


L2(ai  +  a2  +  a3)  +  L3ax 


dL3 

-  =  L2(a3  -  a^)  +  L3(2a3  +  a2  -  a^)  +  L^(a3  +  a2) 

dx 


■ 


. 


- 


. 
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dL^ 

dx 

dL, 


dx 


-L3a2  -  L4(ai  +  a2  -  aj 
L3a2  +  Val  "  a2  '  a3J 


Again  the  performance  criterion  $  is  defined  as: 
5 

$  =  S~  (R.  -  R)2 

i=l 

where  R  -  (  -i - 2 - 2 - i - 5  > 

5 


(2-30) 


(2-31) 


As  in  the  last  example,  the  method  of  steepest  descent  will  be 
utilized  to  ensure  that  $  goes  rapidly  to  its  minimum.  This  requires  that 

2i<  o. 

dx 


But 


2 

T~  “  T(atKi  +  aoKo  +  a.,K0) 
dx  511  22  33 


(2-32) 


where  K, 


(4RxR2  -  4RxR3  +  4R2  -  R2R4  -  R2R5  -  4R2  +  R^ 


+  RqR_  -  5R7  +  5R.R,) 
3  5  4  4  5 


2  2 
K2  -  (-5R  R2  -  R1R3  +  RxR4  +  5R2  -  R2R3  +  R^  +  4R3 


+  4R_RC  -  4R,  +  4R, R,.) 
3  5  4  4  5 


.  -If.  - 


' 


. 


SJJ > 
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■> 
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K3  =  (-6R1R2  -  2R1R3  -  R^  +  4R2  +  2R2R3  -  2R£R4  -  R^ 


+  8R^  +  2R„R.  -  2R_Rt.  +  4R?  -  6R.R,) 
3  34  35  4  4  5 


If  we  let 


al  "  "  K! 


a2  =  “  K2 


(2-33) 


a3  "  -  k3 


then  equation  (2-32)  becomes 


d$  _  2  .  2  2  '  2v 

dx  5  ("al  a2  a3 


It  is  easy  to  see  that  this  is  always  negative. 

Substituting  for  a^,  a2  and  a3  in  equations  (2-29)  and  (2-30) 
using  equation  (2-33)  the  result  in  vector  form  is: 


where 


K  = 


dx 


(2-34) 


dL, 

— D 

dx 


K  L, 

- -D 


(2-35) 


0  K3  +  K2  -  Kx 
0  -(Kx  +  K2  +  K3) 

0  Ki  -  k3 

0  0 

0  0 


K„ 


-K, 


K,  -  Kr 


K, 


-K, 


0 

0 

-  (k2  +  k3) 

K1  +  k2  -  K3 

K2  +  K3  -  Kx 


■ 
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Fig.  2-8  Variation  of  Resistive  Elements  with 
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Fig.  2-9.  Variation  of  Inductive  Elements  with  x. 
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ORIGINAL 

NETWORK 

x  —  0 

NEW 

NETWORK 

x  =  0.002 

R1 

3 

3.414 

R2 

4 

3.572 

r3 

6 

5.013 

R4 

3.5 

4.544 

R5 

7 

6.828 

L1 

1 

1.322 

L2 

3 

2.670 

L3 

5 

3.552 

L4 

8 

9.452 

L5 

6 

6.486 

$ 

11.1 

7.573 

TABLE  2-3  Summary  of  Results 
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414  1.322  5.013  3.552  6.486  6.828 
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Fig.  2-10  Equivalent  Network  at  x  =  0.002 
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The  nonlinear  equations  (2-34)  and  (2-35)  are  solved  as  before  on 


the  digital  computer.  The  solution  is  displayed  in  Appendix  V.  The  graphs 
showing  the  variations  of  the  elements  with  x  are  shown  in  Figs.  2-8  and 
2-9.  The  solution  obtained  for  x  =  0.002  is  not  the  ideal  optimum,  since 
$  can  be  reduced  further.  However,  the  computation  was  stopped  at  this 
value  of  x  since  it  was  felt  that  the  distribution  of  the  resistive  values 
was  good  enough  and  the  values  of  the  other  elements  were  acceptable  for 
purposes  of  this  thesis. 


Table  2-3  shows  the  comparison  between  the  initial  values  and  the 


new  values  of  the  elements*  and  Fig.  2—10  shows  the  new  network. 


The  short  circuit  admittances  for  the  original  and  new  networks 


are : 


Original  Network 


Y(s) 

1-1 


0.311  s2  +  0.586  s  +  0.252 

s3  +  3.684  s2  +  4.181  s  +  1.436 


Y(s)  =  Y (s)  = 
12  21 


2 

0.047  s  +  0.083  s  +  0.028 
s3  +  3.684  s2  +  4.181  s  +  1.436 


Y(s) 

2-2' 


0.107  s2  +  0.276  s  +  0.153 

s3  +  3.684  s2  +  4.181  s  +  1.436 


. 


■ 


.  ,  -  .  :  .•  L  \  r;  »  i 


v 
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Equivalent  network 


/  \  _  0.305  s2  +  0.572  s  +  0.244 

v  /  o  O 

1- 1'  s  +  3.681  s  +  4.175  s  +  1.430 

/  s  v t  \  0.048  S2  +  0.087  s  +  0.031 

(s)  =  Y  (s)  =  —7: - - - 

12  21  s  +  3.681  s  +  4.175  s  +  1.430 

_  0.105  s2  +  0.271  s  +  0.150 

(s)  t  _ 

2- 2'  s  +  3.681  s  +  4.175  s  +  1.430 


Table  2-4  shows  the  percentage  deviation  of  the  coefficients  of 
the  powers  of  s  in  Y(s)  of  the  new  network  from  the  corresponding  coef¬ 
ficients  of  the  original  network. 


%  DEVIATION 

COEFFICIENT 

NUMERATOR  OF 

DENOMINATOR 

OF 

Y  r 

11 

Y12 

Y 

22 

3 

s 

0 

2 

s 

1.93 

-  2.57 

1.87 

0.0815 

1 

s 

2.49 

-  5.32 

1.81 

0.144 

0 

s 

3.175 

-13.23 

1.96 

0.418 

TABLE  2-4 


Percentage  Deviation  in  Coefficients 


“  :  *  '  J  -  ■  Kjl 
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The  large  discrepancy  in  the  case  of  the  coefficients  of  the 

powers  of  s  in  the  transfer  function  seems  to  result  from  the  fact  that 

in  both  the  original  and  equivalent  networks  the  coefficients  of  s  and 
0 

s  are  very  small.  Consequently  even  a  small  error  in  their  values  is 
magnified  when  the  percentage  deviation  is  computed,  since  division  by  a 
small  number  is  involved . 

The  question,  however,  might  be  asked.  Why  not  compare  the  open- 
circuit  parameters?  If  the  open-circuit  parameters  are  considered,  we 
shall  be  dealing  with  at  least  one  less  than  the  number  of  nodes  we  start 
with,  since  1^  or  will  have  to  be  assumed  equal  to  zero.  Hence  in  order 
to  take  into  account  all  the  elements,  1^  and  should  be  nonzero.  One 
way  of  keeping  1^  and  nonzero  is  to  short-circuit  one  of  the  ports.  If 
this  is  done,  then  it  is  more  convenient  to  talk  of  short-circuit  admit¬ 
tances.  This  is  also  according  to  standard  practice. 

2-4  Conclusion 

It  must  be  pointed  out  that  the  equivalent  network  which  is  chosen 
because  it  is  optimal  for  the  spread  in  the  values  of  the  resistances  may 
not  be  desirable  from  other  points  of  view.  For  instance  in  the  two 
examples  considered,  we  may  find  that  in  the  "optimal  network"  the  values  of 
inductances  may  be  too  large.  If  the  designer  is  confronted  with  such  a 
situation,  he  should  perhaps  settle  for  a  "sub-optimal"  network  in  which  the 
distribution  of  resistive  values  is  "good"  and  the  values  of  inductances  are 
also  reasonable.  Another  approach,  a  more  complicated  one,  will  be  to  use 
the  comprehensive  performance  criterion  as  given  by  equation  (2-1) .  By  a 


. 


. 


; 


. 


- 
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■ 
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judicious  choice  of  values  of  the  weighting  factors,  the  relative  importance 
of  one  or  more  kinds  of  elements  can  be  controlled. 


Chapter  3 


NETWORK  OPTIMIZATION— II  (Sensitivity) 


3-1  Introduction ;  Recently  a  number  of  papers  [2,  9,  16]  has  been 
published  on  the  application  of  the  concept  of  continuous  equivalence  to 
minimize  the  sensitivity  of  network  functions  to  changes  in  the  element 
values.  In  this  approach  the  problem  of  minimization  of  the  sensitivity  of 
a  network  consists  of  scanning  through  a  set  of  networks  which  are  all 
equivalent ,  and  then  picking  one  with  minimum  sensitivity.  A  summary  of 
Schoeffler’s  approach  to  this  problem  is  given  in  Appendix  III.  It  is  felt 
that  a  simpler  approach  can  be  used  in  most  cases  except  perhaps  when  the 
number  of  elements  is  unduly  large.  This  procedure  will  be  described  in 
section  two  of  this  chapter.  Before  we  do  that  however,  a  brief  discussion 
on  the  classical  definition  of  sensitivity  will  be  presented. 


Let  T  be  a  network  function,  e.g.,  driving  point  impedance.  If 


we  represent  the  N  elements  in  the  network  by  e^,  e^t 


eK* 


then  T  can  be  expressed  as 


T  =  T (s ,  e±,  e2. 


where  s  is  Laplace  transform  variable  or  complex  frequency  and 


e  is  any  element  (L,  R  or  C) . 
K 


Using  the  classical  definition  of  sensitivity,  the  sensitivity 

th 

S  ,  of  T  to  a  change  in  the  K  element  of  the  network  is  given  by 
K 


SK  3 in  eK  T  3eK 


3 In  T  eK  3T 


1 1  I  Hllllll  I II 
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3-2  Minimization  of  Sensitivity 

Let  us  consider  the  one-port  network  in  Fig.  3-1  which  has 

+  N  +  N  elements.  (N  ,  N  and  N_  are  the  number  of  resistors, 
k  L  C  R  L  C 


9 


Fig.  3-1  Passive  R-L-C  Network 

inductors  and  capacitors  respectively.)  Let  us  assume  that  we  are  required 
to  find  its  equivalent  in  which  the  sensitivity  is  minimum.  Since  the 
sensitivity  function  S  is  complex,  a  convenient  and  meaningful  perform- 
ance  criterion  as  suggested  by  Schoeffler  [16]  can  be  expressed  as 


«  . .  ,  ■  l  :  a 


• 
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$ 


WNC 

II 


K=1 


(3-2) 


2  NL  2  nc  2 


*  L 

+  y~ 

s 

*■  C 

+  Y 

s  , 

Ri 

2_ 

j=i 

Lj 

4=i 

Ck 

(3-3) 


th 


where  S  is  the  sensitivity  of  T  to  a  change  in  the  i  resistive 
element.  S  and  are  similarly  defined. 


Our  purpose  then  is  to  find  an  equivalent  network  which  has 
minimum  value  of  $. 

As  in  the  previous  examples,  the  optimization  problem  is  solved 
in  two  stages.  First  the  generating  equations  are  obtained  and  then  $ 
is  minimized  using  the  generating  equations  as  constraints. 

The  procedure  for  deriving  the  generating  equations  has  been 
discussed  earlier  in  Chapter  1  (equations  1-22,  1-23,  1-24).  We  can  there¬ 
fore  write  the  generating  equations  for  the  network  represented  by  Fig.  3-1 
as 


dR. 

— D 
dx 


dx 


M  L, 
- b 


(3-4) 


dx 


!  ' 


' 
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The  order  of  M  will  depend  on  the  actual  network  configuration,  number  of 
branches,  etc.  The  elements  of  M  are  generally  linear  combinations  of 
the  nonzero  elements  of  A^. 

For  simplicity,  let  N  =  N  =  N  =  N  so  that  equation  (3-2) 

R  L  C 

becomes 


4> 


(3-5) 


As  before,  the  method  of  steepest  descent  will  be  used  to  ensure  that  $ 

d$ 

goes  rapidly  to  its  minimum.  This  requires  that  <  0 


d$ 

dx 


de 


K 


dx 


<  0 


(3-6) 


where  e^  is  any  element  (R,  L,  or  C) . 

deK 

If  we  substitute  for  -j^—  ,  using  equation  (3-4),  in  equation  (3-6), 
we  shall  get  an  expression  of  the  form 


d£ 

dx 


2_ 


)  <  0 


(3-7) 


where  n 


is  the  number  of  nonzero  elements  of  A^  which  form  the  matrix  M. 
By  defining  a±  =  -f±(e1,  e2>  .  3^)  .  (3-8) 


and  substituting  in  (3-7) ,  we  get 


n 


•  I 


-  V 


i=l 


d$ 

dx 


(3-9) 


.  ...  .  ’  •• 
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It  is  easy  to  see  that  is  always  negative. 

If  we  then  substitute  for  a^’s  as  defined  by  equation  (3-8)  into 
equation  (3-4) ,  we  can  solve  the  resulting  nonlinear  equation  on  a  digital 
computer.  This  method  of  solution  is  more  straightforward  than  Schoeffler's 
approach.  It  also  has  the  following  advantages:  (1)  Only  one  set  of  differ¬ 
ential  equations  will  be  solved  instead  of  two  as  Schoeffler’s  technique 
requires.  (2)  Since  the  solution  is  self-adjusting,  the  problem  of  making 
an  initial  and  consequent  choices  of  a^'s,  every  time  the  generating  equa¬ 
tions  are  solved,  is  eliminated.  (3)  The  problem  of  putting  bounds  on  a^’s 
is  also  eliminated. 

In  the  next  section,  an  example  which  illustrates  the  optimization 
procedure  discussed  above  is  presented. 

3-3  Example  3-1 

Consider  the  minimization  of  the  sensitivity  of  the  network  shown 
in  Fig.  3-2.  Since  the  structure  of  the  network  is  the  same  as  in  Example 
2-1,  the  element  generating  equations  are  the  same  as  equations  (2-8),  and 
(2-10) .  The  equations  are  repeated  below  for  convenience. 

dR, 

—  =  -  R2(a2  +  ai>  +  R3al 

dx 

dR 

-  =  R2^a2  “  al^  ~  R3ai  .  (3-28) 

dx 

dR„ 

— -  =  R2(a£  +  a^  +  R3(2a2  +  a^ 

dx 


.  •  • 

* 

•  ?  u 

abnuoc 
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R;l=1  h1=l  R3=l 


Fig.  3-2  Given  R-C  Passive  Network 


-  h2(a2  +  a1)  +  h3a1 


dh2 

=  h2(a2  “  a^)  “  ^3ai 

dx 

dh3 

-  =  h2  (a2  +  a1)  +  h3(2a2  +  a^ 

dx 


where  as  before 


h 


i 


(3-29) 


v 


. 
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We  note  in  Fig.  3-2  that  there  is  no  capacitor  in  the  third  branch, 
consequently  h^  =  0. 

In  order  to  keep  the  cost  down  it  is  our  intention  to  keep  the 
number  of  elements  the  same  as  in  Fig.  3-2.  In  other  words,  we  shall  keep  h^ 
at  zero  in  the  equivalent  network  also. 

Before  we  take  up  the  minimization  of  sensitivity,  we  will  point 
out  the  constraints  which  have  to  be  imposed  on  the  generating  equations  in 
case  an  element  is  zero  to  begin  with  and  we  wish  to  keep  it  at  zero  or  one 
of  the  elements  becomes  zero  in  one  of  the  equivalent  networks  and  we  do  not 
want  it  to  become  negative.  The  constraints  consist  of  equating  the  particu¬ 
lar  element  and  its  derivative  with  respect  to  x  to  zero  in  the  generating 

dh^ 

equations.  Thus  in  equations  (3-29)  if  h^  =  0  and  =0  we  get 

dhi 

-  =  -  h2(a2  +  ax) 

dx 

dh9 

—  =  h2(a2  -  ax)  .  (3-30) 

dx 

h2(a2  +  a1)  =  0 


Since  h2  ^  0 


a 


2 


(3-31) 


Equation  (3-31)  gives  the  additional  constraint  which  should  be  imposed  on 
equations  (3-28)  and  (3-29)  if  equivalent  network  is  to  result  in  which 


0. 


. 


' 
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get 


Substituting  equation  (3-31)  into  equations  (3-28)  and  (3-29)  we 


dR., 


dx 


R3  al 


dR, 


dx 


2^2  ai  ”  ^3  a]_ 


dR, 


dx 


-  R3  a1 


dx 

dh, 

dx 


=  0 


-  2h2 


(3-32) 


Now  we  proceed  with  the  minimization  of  sensitivity.  If  we  now  define  the 
sensitivity  criterion  $  as  in  equation  (3-1)  our  problem  is  to  minimize 
$  using  equations  (3-32)  as  constraints. 

For  the  example  under  consideration,  $  is  given  by 


$ 


S 


i 


3Z  . 

3e .  Z 

l 


(3-33) 


where 


B 


I  i  |  8. 

. 
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th 

e_^  is  the  i  element  R  or  h 


an<3  Z  is  the  driving  point  impedance, 


A  calculation  of  $  for  the  network  discloses  that  the  series 
arm  or  branch  one  contributes  89.3 %  of  the  value  of  $.  This  confirms 
Milkurski's  [11]  claim  that  elements  in  series  are  more  effective  in  their 
change  than  those  in  parallel.  Thus,  the  effect  of  an  element  variation  on 
a  network  function  depends  a  great  deal  on  the  position  of  the  element  in 
the  network. 

In  view  of  this,  a  simplified  version  of  equation  (3-33)  can  be 
used  such  that  only  the  series  arm  is  taken  into  consideration.  The  new 
performance  criterion  is  designated  and  expressed  as: 


(3-34) 


By  finding  an  equivalent  network  with  a  reduced  value  of  it 

is  expected  that  $  will  be  reduced.  However,  due  to  the  nonlinear  nature 
of  $,  it  will  not  change  linearly  with 

The  remaining  computations  will  be  carried  out  at  a  constant  value 
of  frequency,  namely,  co  =  1  radian  per  second.  (In  other  words,  we  have 
let  s  =  jo)  =  jl  in  equation  (3-34).) 

Equation  (3-34)  can  be  written  as 

(R2  +  h2)  {(R2  +  R3)2  +  hj} 

$1  - 

{[h.h2  -  (RjR,  +  R.R3  +  R2R3)]2  +  (hlR2  +  hlR3  +  h2Rl  +  h2R3)2} 


(3-35) 


* 


' 
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For  simplification  purposes,  the  following  definitions  are  made 


K11 

= 

(r2  -1 

■  R3)2  +  h2 

K12 

= 

R2  + 
1 

h2 

1 

K13 

= 

hlh2 

-(R^R2  +  R^R3  +  R2R3^ 

K14 

- 

hlR2 

+  h^R3  4-  ^2Ri  ^2R3 

K15 

= 

2  2 

K13  +  K14 

K16 

= 

K11  K12 

K17 

= 

2Kn 

K15  R1 

K18 

s 

2K16 

{Kuh2  -  K13(R2  +  R3)} 

K19 

= 

2K11 

K15  hl 

K20 

= 

2K16 

{K13h2  +  K14(R2  +  R3)} 

K21 

= 

2K12 

k15(r2  +  r3) 

K22 

= 

2K16 

^K14hl  "  K13^R1  +  R3^ 

K23 

= 

2K16 

{K14(hl  +  h2)  "  K13(R1  +  R2)} 

K24 

= 

2K12 

K15  h2 

K25 

= 

2K16 

{K  h.  +  K  ,(1^  +  R3)} 

. 


. 
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K26 

-  r3(k1?  + 

K22  ' 

K27 

=  2R2(K22  - 

K21 

K29 

=  2h2(K25  - 

K24 

K30 

K26  +  K27 

+  K, 

ir 


(3-36) 


Applying  the  method  of  steepest  descent  on  to  arrive  at  its 


minimum  requires 


d<&. 


<  0 


.  (3-37) 


dx 


Substituting  equation  (3-35)  in  (3-37)  and  using  the  definitions  in  (3-36) 
the  result  is: 

K, 


a,  •  -f  <  0 

KL 


(3-38) 


One  way  of  satisfying  inequality  (3-38)  is  to  define  a^  as: 


-K 


al  = 


30 


(3-39) 


K 


15 


Although  is  obviously  positive,  it  is  retained  in 

equation  (3-39)  because  (as  it  was  discovered  during  numerical  computa- 

2 

tions)  K30  may  be  too  large  for  the  computer  to  handle.  Retaining 


' 
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- 


_ 
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enables  us  to  control  the  value  of  a^  since  like  K^q  depends  on 

the  values  of  the  elements  of  the  network. 

Substituting  for  a^  in  equation  (3-32)  the  final  generating 
equations  are : 


dR„ 


dx 


K 


=  -  R 


30 


3  2 
K15 


dR, 


dx 


2K 


=  R„ 


K 


30  R 
2  +  r3 
15 


K 


30 


K 


15 


dR, 


dx 


=  R, 


K 


30 


K 


15 


dh. 


dx 


=  0 


dh. 


dx 


=  h 


2K 


30 


2  K2 
K15 


(3-40) 


A  part  of  the  solution  of  equation  (3-40)  is  displayed  in 
Appendix  VI  and  the  curves  representing  the  variations  of  the  elements  with 
x  are  shown  in  Fig.  3-3. 

Assuming  that  we  are  satisfied  with  the  value  of  $  at  x  =  0.177, 
the  network  that  results  is  shown  in  Fig.  3-4. 
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0.781  1.0 


1.219 


Fig.  3-4  Equivalent  Network  at  x  =  0.177 
R  is  in  ohms  and  C  is  in  inverse  farad. 


Table  3-1  shows  the  comparison  of  the  different  network  elements 
and  their  respective  $  at  to  =  1. 

The  driving  point  impedances  of  the  original  and  the  new  networks 


are : 


Original  Network  (Fig.  3-2) 

v  S2  +  1.333  s  +  0.333 
Z  (s)  =  - 2 - 

0.667  s  +  0.333  s 

New  Network  (Fib.  (3-3) 

s2  +  1.325  s  +  0.337 

0.625  s2  +  0.337  s 


Z(s) 


' 


' 

' 


. 
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ORIGINAL 

EQUIVALENT 

NETWORK 

NETWORK 

o 

ii 

X 

x  =  0.177 

R1 

1 

0.781 

r2 

1 

1.753 

R3 

1 

1.219 

hl 

1 

1.0 

h2 

1 

1.534 

h3 

0 

0 

$ 

0.56 

0.51 

TABLE 

3-1  Summary  of 

Results 
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Table  3-2  shows  the  percentage  deviation  of  the  coefficients  of 
the  powers  of  s  in  the  Z(s)  of  the  new  network  from  the  corresponding 
coefficients  of  the  original  network. 


COEFFICIENTS 

%  DEVIATION 

OF 

NUMERATOR 

DENOMINATOR 

2 

s 

0 

2.1 

1 

s 

0.6 

-  1.333 

0 

s 

-  1.333 

— 

TABLE  3-2  Percentage  Deviation  in  Coefficients 

It  should  be  mentioned  at  this  point  that  the  new  network  repre¬ 
sented  at  x  =  0.177  is  by  no  means  the  minimum  sensitivity  network  since 
the  value  of  $  can  be  reduced  further.  However,  the  solution  was  stopped 
at  x  =  0.177  since  the  main  objective  of  this  discussion  is  to  establish 
the  usefulness  of  the  procedure.  As  stated  earlier,  the  choice  of  a  net¬ 
work  depends  on  several  factors  and  the  designer  will  have  to  make  the 
decision.  Some  of  these  factors  are  discussed  in  Chapter  4. 
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3-4  Comments 

Although  the  minimization  of  sensitivity  was  carried  out  with 
respect  to  the  elements  in  the  series  branch,  this  procedure  can  equally 
well  be  used  to  minimize  the  sensitivity  with  respect  to  any  other  elements 
or  groups  of  elements  in  the  network. 

Schoeffler  [16]  and  Leed  [9]  have  considered  examples  which  lead 
to  the  conclusion  that  the  equivalent  network  with  minimum  sensitivity  has 
more  elements  than  the  original  network.  If  their  procedure  is  applied  to 
the  example  discussed  in  this  chapter,  the  minimum  sensitivity  network  may 
turn  out  to  be  one  with  at  least  six  elements.  The  minimum  value  of  sensi¬ 
tivity  of  the  six-element  network  may  be  less  than  the  minimum  value  for  the 
five-element  network.  Even  if  this  is  so,  it  must  be  remembered  that  the 
advantage  gained  by  reducing  sensitivity  may  be  offset  by  the  increase  in 
cost  due  to  the  additional  elements.  A  judgement  will  have  to  be  made  in 
each  individual  case. 


. 


- 
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Chapter  4 


SUMMARY  AND  CONCLUSION 


Some  Limitations  of  Continuous  Equivalence  Theory 

Although  Schoeffler  [15]  published  his  theory  of  continuously 
equivalent  networks  in  1964,  only  during  the  past  couple  of  years  have 
reports  on  some  of  the  limitations  of  the  theory  appeared  in  the  literature. 
In  1966,  Newcomb  [12]  produced  an  example  of  two  networks  with  the  same  con¬ 
figuration,  the  same  number  of  elements  and  the  same  driving  point  impedance. 
He  showed  that  neither  of  these  networks  could  be  obtained  from  the  other  by 
Schoeffler's  method.  Thus,  he  concluded  that  not  all  the  equivalent  networks 
can  be  found  by  Schoeffler's  method.  Newcomb  arrived  at  this  conclusion  by 
solving  the  generating  equations. 

By  way  of  amplifying  Newcomb's  statement,  it  must  be  pointed  out 
that  other  equivalent  networks  of  either  of  the  two  networks  considered  by 
Newcomb  can  be  generated  by  Schoeffler's  method  provided  we  are  allowed  to 
increase  the  number  of  elements.  Another  point  which  may  be  worth  mention¬ 
ing  is  that  it  is  not  always  necessary  to  solve  the  generating  equations,  as 
Newcomb  did,  to  conclude  whether  Schoeffler's  theory  fails  or  not  in  a 
particular  case. 

In  order  to  illustrate  this  point,  consider  the  following  example: 

Example  4-1 

Suppose  we  are  given  the  network  in  Fig.  4-1  (this  is  the  same 
network  as  in  example  2-1  except  that  R£  =  0  and  h^  =  0)  and  we  are  asked 
to  find  an  equivalent  network  with  the  same  configuration  and  same  number 
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of  elements. 


R1  L1  C1  R3  L3 


» — vw — w— |H 

— vw1 — 'Wp — 

5  L0 

3  2 

-  Cn 

-v  2 

Fig.  4-1  One 

-Port  R-L-C. 

Dimensions  are  in  ohms,  henry  and  daraf . 


We  can  use  the  generating  equations  obtained  for  example  2-1  but  must  allow 
for  the  fact  that  =  0  and  h^  =  0.  The  equations  are: 


dx 


-  R2(a2  +  a2)  + 


R2 (a2  ”  al^  —  ^3ai 

R2  (a2  +  al^  +  R3^a2  +  al^ 


dx 


■ 


i 
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dx 

dL, 

dx 

dL, 

dx 

dx 

dh, 

dx 

dh, 

dx 


-  L2^a2  +  al^  +  L3al 


L2^a2  al^  “  L3al 


^2^a2  +  al^  +  ^3^a2  +  al^ 


-  h2(a2  +  a1>  +  h3a1 


h2(a2  *  al)  '  h3al 


h2(a2  +  ax)  +  h3(2a2  +  a.^ 


(4-1) 


Since  R2  =  0  and  h3  =  0,  we  impose  the  constraints 


dR, 


dx 


dh. 


dx 


=  0;  R2  -  0 


=  0;  h3  =  0 


(4-2) 


on  equation  (4-1) . 

The  generating  equations  become 


dR-j 

dx 

dR, 

dx 


=  0 


=  0 


- 


' 


. 
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dL: 

dx 

dL, 

dx 

dL, 

dx 

dh^ 

dx 

dh, 

i 

dx 


=  0 


=  0 


=  0 


=  0 


=  0 


(4-3) 


Since  equations  (4-3)  imply  that  all  the  elements  are  constant,  it 
can  be  concluded  that  Schoeffler's  method  will  not  yield  a  network  equivalent 
to  Fig.  4-1,  if  the  same  configuration  and  number  of  elements  is  to  be  main¬ 
tained.  It  may,  however,  be  possible  to  find  an  equivalent  network  if  the 
number  of  elements  can  be  increased  since  this  will  mean  the  deletion  of  the 
constraints  given  by  equation  (4-2) . 

Another  situation  where  Schoeffler's  method  fails  is  related  to 
multiport  networks.  In  view  of  the  fact  that  Schoeffler's  method  is  an 
extension  of  Howitt's  work,  one  of  the  conditions  which  must  be  satisfied  for 
Schoeffler's  method  to  be  valid  is  that  the  number  of  loops  must  be  greater 

than  the  number  of  ports.  There  is  no  provision  in  Schoeffler's  method  to 
alter  the  number  of  loops  in  the  course  of  the  solution.  For  instance, 
Schoeffler's  method  cannot  be  used  to  find  an  equivalent  of  the  two-port 


-  65  - 


network  shown  in  Fig.  4-2  because  the  network  has  only  two  loops. 


Lastly,  the  network  shown  in  Fig.  4-3  is  another  example  in  which 
Schoeffler's  theory  does  not  yield  an  equivalent.  The  reason  for  the  non¬ 
equivalence  for  this  case  is  still  under  investigation.  Analytically,  the 
theory  is  supposed  to  yield  an  equivalent,  but  computer  solution  has  revealed 
that  the  deviation  of  the  driving  point  impedance  of  the  equivalent  network 
from  that  of  the  original  network  increases  as  x  increases. 


Fig.  4-3  One-Port  R-C  Network. 
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-• 


'  ■  if. 


-  66  - 


4_2  Comments  on  the  Choice  of  Optimal  Networks 

In  all  the  examples  discussed  in  this  thesis,  caution  has  been 
exercised  in  referring  to  any  particular  solution  as  the  optimal  solution. 

The  main  reason  is  that  a  network  which  is  optimal  with  respect  to  one 
performance  criterion  may  prove  to  be  undesirable  from  other  points  of  view. 
The  choice  of  a  truly  optimal  network  depends  largely  on  practical  consider¬ 
ations.  Some  of  these  considerations  are:  small  values  of  capacitances  are 
undesirable  since  they  are  difficult  to  measure  accurately;  large  values  of 
inductances  are  affected  by  shunt  capacitances  between  their  windings.  It 
is  therefore  left  to  the  network  designer  to  exercise  his  discretion  based 
on  practical  experience  to  choose  what  he  considers  to  be  the  optimal 
network. 

4-3  Areas  for  Further  Study 

As  pointed  out  earlier,  Schoeffler’s  method  may  not  yield  all  the 
equivalents  of  a  given  network  or  may  yield  none  at  all.  At  the  present 
time,  in  order  to  determine  whether  Schoeffler's  method  will  work  or  not, 
one  has  to  obtain  the  generating  equations  which  in  itself  may  be  a  time- 
consuming  task  especially  if  the  network  has  several  elements.  Therefore  it 
would  appear  that  one  of  the  areas  for  further  research  should  deal  with  the 
development  of  simple  procedures  or  tests  to  determine  whether  the  equiva¬ 
lent  of  a  given  network  can  be  found  by  Schoeffler’s  method. 

If  the  limitations  of  the  method  can  be  eliminated,  Schoeffler's 
continuous  equivalence  theory  may  prove  to  be  one  of  the  best  tools  available 
to  network  designers  for  optimization  of  networks.  Hence,  further  research 
is  worth  while  to  eliminate,  if  possible,  at  least  some  of  the  limitations. 
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APPENDIX  I 


Theorem: 
and  Z^(s)  have  the 


Two  (n  4-  m)  x  (n  +  m)  loop- impedance  matrices 
same  n  x  n  terminal-impedance  matrix  Z(s)  if 


Z  (s) 
a 


Z  =A  Z  A 


where 


A  = 


u  1  ot, 

-n  |  “12 

o  ;  oc 

t  I 

n  m 


}  n 


} 


m 


A  =  • 


u  I  0 

— n  .  — 


—21  I  —22 


}  n 

}  m 


n 


m 


and  3\j  is  assumed  nonsingular. 


(Al-1) 


(Al-2) 


.  (A1  -  3) 


Proof : 


final  m 


Consider  the  (n  +  m)-port  described  by  Z^. 
ports  yields  the  description 


Short-circuiting  the 


"h 

}  n 

-  4 

M|  1 

1 _ L 

=  Jl  ZaA 

~h 

0 

}  m 

--2- 

-h. 

(Al-4) 


where  the  voltage  and  current  vectors  are  partitioned  as  the  ports. 


■ 


■ 
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Premultiplying  equation  (Al-4)  by  J\j  1  we  get 


V, 

I, 

-1 

=  z 

-1 

~a 

0 

T3 

Ml 

_ 1 

where 


-2 


-211!  +  -22I2 


Since  the  variables  and  1^  for  all  1^  are  the  same  for  the 

two  circuits,  they  must  be  described  by  the  same  terminal  impedance  matrix. 


?1 


-ZI-i 


Q.E.D. 


The  theorem  and  its  proof  is  taken  from  Newcomb  [13] ,  but  the 
notations  have  been  changed  to  conform  with  the  notations  in  this  thesis. 
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APPENDIX  II 


A  summary  of  Schoeffler’s  [15]  development  of  continuously 
equivalent  networks  is  presented  in  this  appendix.  Some  of  the  notations 
have  been  changed  for  clarity. 

In  order  to  formulate  the  problems  of  continuously  equivalent 

networks,  consider  the  transformation  matrix  which  transforms  the  branch 

admittance  matrix  Y,  into  Y,  +  AY,  where  AY,  is  incremental  in  size. 

— b  — b  --b  — b 


Let  T  = 

u  +  AT 

.  (A2-1) 

where  u  is  b  x  b 

unit  matrix  and 

b 

is  the  number  of  branches. 

The  application  of 

T  to  Y,  results 
—  — b 

in 

Y^  given  by 

ii 

T?  Y,  T 
—  — b  — 

.  (A2-2) 

= 

Y,  +  Y^AT  +  ATfc 
— b  — b  — 

4 

+  ATC  Y,  AT  .  (A2-3) 

—  — b 

define 

n 

T^Ax 

.  (A2-4) 

where  Ax 

is  a  scalar. 

- 

Then  T  and 


Y,  become 
— b 


. 


. 


■ 
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Writing  T  =  T- (x)  ,  Y,  =  Y,(x)  and  Y7  =  Yu(x+Ax) 

— ±  — X  — D  — D  — D  — D 

equation  (A2-6)  becomes 


Y^  (x  +  Ax)  -  Y^x) 
Ax 


Y  T  +  T  Y  + 
-1-1  -b 


0(Ax2) 


Ax 


(A2-7) 


Passing  to  the  limit  yields  the  matrix  differential  equation 


dY 

-  =  Yb  Tx  +  tJ  Yb  .  (A2-8) 

dx 

The  matrix  T-^(x)  must  be  selected  so  that  equivalent  networks 
result.  To  derive  the  equivalence  restrictions  on  T-^,  let  us  start  with 
the  general  restriction 

Ta11  =  afcA  .  (A2-9) 

where  a  is  the  usual  n  x  b  cut  set  matrix  of  the  network,  A  is  Cauer 
transformation  matrix,  and  n  is  the  number  of  independent  node  pairs. 

Writing 


T  =  u  +  TxAx  .  (A2-10) 

and 

A  =  u  +  BAx  .  (A2-11) 


equation  (A2-9)  becomes 


P, 

. 

■ 


. 

■ 


J 


1 


. 
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t  t  t  .  t 

a  +  T^a  Ax  =  a  +  a  BAx 


or  Tafc  =  atB 


(A2-12) 

(A2-13) 


where  B^  is  of  the  form  (considering  the  preservation  of  the  driving  point 
and  transfer  impedances  at  terminal  pairs  1  and  2) 


B  = 


0 

0 

0 

0 

0 

0 

0 

0 

X 

X 

X 

X 

X 

X 

X 

X 

(A2-14) 


For  the  most  general  branch  in  which  R,  L  and  C  are  connected 
in  parallel, 


Y,  =  G  +  sC  +  -  L  1 
— b  -  s  - 


(A2-15) 


where  £  is  the  branch  conductance  matrix 
C  is  the  branch  capacitance  matrix 
and  L  ^  is  the  branch  inverse  inductance  matrix, 


Thus,  the  transformation  T  can  be  applied  to  each  element  kind  separately 
provided  the  same  Cauer  transformation  matrix  A  results  for  each  element 
kind.  That  means  the  branch  admittance  matrices  may  be  transformed  as 


£'  -  Soifiloi 


£  ■  Jci£-Tci 


<t1> '  “  Ili  l-'X  Sli 


(A2-16) 

(A2-17) 

(A2-18) 


-G1 "  Tci^  =  "Ll”'  =  a-S- 


prov ided 
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APPENDIX  III 


Schoeffler  [16]  presents  an  approach  to  the  synthesis  of  minimum 
sensitivity  networks.  Below  is  a  summary  of  the  relevant  equations  which 
lead  to  the  minimum  sensitivity  network. 

The  generating  equations  are  generally  of  the  form 


d£ 

dx 


=  MG 


d£ 

dx 


=  MC 


(A3-1) 


dL 

dx 


—  =  ML 


where  G  is  the  branch  conductance  matrix 
C  is  the  branch  capacitance  matrix 
and  L  is  the  branch  inductance  matrix. 

Define  the  sensitivity  of,  say,  the  transfer  function  T  to  a 
change  in  the  k^  element  e,  by 


q  =  9T_ 

k  T  *  3e, 


(A3-2) 


Also  define 


sk  =  \  ek 


(A3-3) 


After  some  mathematical  manipulations,  Schoeffler  arrives  at  the  equations 
for  all  the  sensitivities  in  matrix  form  given  by: 
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f2c 

dx 


dx 


*  -  M 


M  Qg 


-«  9c 


(A3-4) 


dx 


=  -  M  Qt 


where 


-G  *-qGl’  qG2 


qGb^ 


Qc  anc^  Ql  can  be  similarly  defined,  and  t  indicates  the  transpose, 
Define  the  performance  criterion  $  by: 


$ 


-  T  H  [SGk  SSk  +  SCk  SSk  +  SLk  SLk] 


(A3-5) 


where  *  indicates  conjugate. 


£  ‘  -9g  [5g  ac  +  gt  V  -QG  « 

-  4  [p l  MC  +  MD*]  Q*  +  Cc  Dqc  D*c  MC  (A3-6) 

-  2 *  tsj  MC  +  ffl*]  Q*  +  LC  DQL  D*L  ML 

where  the  Q  ,  Q  ,  Q  ,  G,  C  and  L  are  column  matrices  whose  elements  are 
C?  L*  L< 

the  various  and  the  values  of  the  elements  of  the  network  and  where 

the  D  matrices  are  diagonal  matrices  defined  by: 


Dg  =  Diag  [gx,  g2  . gb] 

?QG  =  Diag  ^qGl*  qG2  .  qGb^ 


■ 
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For  steepest  descent  to  the  minimum,  —  should  be  as  negative  as 
possible.  Thus,  the  free  elements  b^  (|b^|  ^;1)  which  appear  in  M  should 
be  chosen  to  make  this  so. 

d$ 

"The  optimization  of  the  network  then  proceeds  by  determining  — , 
choosing  the  b’s,  calculating  the  transformation  matrix  M,  transforming  to 
another  equivalent  network  by  integrating  equations  (A3-1)  and  (A3-4) ,  choos¬ 
ing  another  set  of  b's,  etc., 


until  the  minimum  is  obtained. 
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